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Abstract 



We describe how the non-local regularization can be implemented in the cal- 
culation of solitons in the Nambu Jona-Lasinio model as well as in the equivalent 
linear c-model. We investigate different forms of regulators and show that the 

!*""■" . 3-momentum cut-off leads to serious conceptual difficulties. 
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1 Motivation 

This work was done together with Georges Ripka and Wojciech Broniowski. 
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D . Solitons corresponding to baryons have been found in several chiral quark models. 

Many of these solutions turn out to be unstable against collapse unless additional con- 
straints are introduced in the model. The well known examples are the linear NJL model 
with proper time regularization |j], [2j and the linear a-model with sea quarks [|3], ||. Even 
in the linear a-model with only valence quarks the energy of the soliton becomes too low 
for any choice of model parameters if one goes beyond the mean field approximation. In 
all these models the instability occurs because it is energetically favorable for the chiral 
field to acquire arbitrary (or very) high gradients. This suggests that cutting off high 
momenta in the interaction may prevent the collapse and stabilize the soliton. A simple 
sharp cut-off does not yield a stable solution while a smooth behavior of the regulator 
(usually interpreted as a fc-dependent quark mass) can indeed lead to solitons which are 
stable against the decay into free quarks as well as against collapse. Such a regularization 
has a physical justification in QCD calculations of the quark propagation in an instanton 
liquid which predict a non-local effective interaction between quarks with a 4-momentum 
cut-off A ~ 600 MeV §. 

Further physical implications of the non-local regularization are discussed in the con- 
tributions to this workshop by George Ripka and Wojciech Broniowski ||. 
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2 The NJL model with non-local regulators 

The non-local regularization of the quark-quark interaction can be implemented in the 
NJL type models by replacing the contact term (q(x)T a q(x)) , T a = (l,i7 5 r a ) by a non- 
local form. Usually one introduces a regulator r diagonal in 4-momentum space such 
that q(x) — > J d^y(x\r\y)q(y). The QCD derivation of the quark propagation in a dilute 
instanton gas predicts the following functional dependence for r(k 2 ) = (k'\r\k)5(k — k') 

§: 

r = f{z) = -z± (Jo (z) K (z) - h (z) Ki{z)) , z= ^ , (1) 

where p is the instanton size of the order (600 MeV)" 1 . As we shall see in the following 
it is necessary to analytically continue the regulator to negative k 2 in order to be able to 
treat the valence orbit. This is not possible with the form ([!]) since it has a cut along the 
negative real axis starting at k 2 = 0. We use instead a Gaussian shape of the regulator: 

r(k 2 )=e-^i (2) 

or a "monopole" shape: 

r(k 2 ) = -V > (3) 

1 "I" A 2 

which has the proper behavior for large k 2 where one gluon exchange dominates. 

The expression for the energy of the soliton and the self-consistency equations can be 
derived from the bosonized Euclidean action 

/ = -Tr log (-ifi^ + m + r(S + i l5 P a r a ) r) + ^ J d 4 x (S 2 + P 2 ) , (4) 

where S and P a are the chiral fields and are the dynamical variables of the system. 

The main difficulty is the presence of time in the regulator. In order to evaluate the 
trace in (^) it is convenient to introduce energy dependent basis states, which are solutions 
of the Dirac equation: 

h ^ 2 )\(lju) =e j (u 2 )\q ju ) (5) 

with 

h (u 2 ) = -ia • V + (3r(v 2 - V 2 ) (S(r) + ij 5 T a P a (r)) r{v 2 - V 2 ) + (3m . (6) 

From (f|) the following expression for a stationary configuration can be derived [0: 



c1c j 



E= Tn L •*> E -^ + ^!^ (S 2 + Z) - v-uum . (7) 

Note that when no regulator (i.e. r = 1) or a time-independent regulator is used, the 
energies Cj are independent of v and the integration can be carried out using the Cauchy 
theorem. Closing the contour from below yields the well known expression for the energy 
of the Dirac sea: E se3, = J2e <o e j- (Note that the energies of occupied orbits lie on the 
negative imaginary z/-axis.) 



When the soliton describes a baryon, the energy of three valence quarks is added to 
the energy of the Dirac sea.[] The same result can be formally obtained by deforming the 
contour in (|7|) in such a way as to encircle the valence orbit (for detailed discussion on 
this point see Wojciech Broniowski contribution to this workshop). Such a prescription 
gives the expected result provided the orbits do not depend on v. However, when the 
regulator depends on time (or z/), this may not lead to the correct result since the regulator 
generates additional poles scattered in the whole complex z/-plane. It may still work well 
for an isolated pole on the positive imaginary axis close to as is the case of the + 
orbit in the soliton with the hedgehog form of the background chiral field ||. This pole 
can then be treated separately, yielding the valence contribution to the soliton energy 
E val = 3e va i, where the energy of the valence orbit is determined from 



iz/ + e +(z/ ) 2 =0. 

val 

The soliton energy can now be written as: 

rp pval _i_ psea _i_ pmeson 



(9) 



The sea contribution is 
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e> 2 )(e> 2 ) - 2v 2 b d (y 2 )) e^VXe.V) - 1v 2 b]{v 2 )) 



v 2 + ej(u 2 ) 2 



v 2 + e){y 2 ) 2 



(10) 



with bj(u 2 ) = dej(v 2 ) / 'dv 2 and is evaluated by direct numerical integration along the 
real z/-axis. The term E meson is given by the last integral in (0) (with the integrand 
S 2 + P 2 -M 2 ). 

The above prescription is further supported by the fact that it gives an exact result 
for the baryon number, which can be expressed as @: 
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dej(u) 
dv 



iz/ + ej(u) 
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The self-consistent equations derived from @ take the form (the hedgehog ansatz, 
Pa{f) — r a P{r), for the pion field is assumed): 
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where <ij v (r) = r((u 2 — V 2 )/A 2 )g Ji/ (r) and res^ = 1 — i-^7 2 - is the residue of the valence 
pole. 

A necessary condition for a stable soliton configuration is that the energy (|7|) is lower 
than the energy of three free quarks. When the regulator depends on time, the free quark 



Only if it is positive, otherwise it is already contained in the above sum. 





-1 




1 


1 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 1 | 1 1 1 1 


i i i i i i i i i— 











[ 


~ 


1800 








1 


i 


1600 








1 








= 


1400 








J ; ,•"'*' ^ so1 


-E 










/. y* 


E 


1200 








* ■ — 


•i 










>^7 


z 


1000 








/ 


-E 








y 




- 






y 






= 




c. / 


\ 






z 


800 






\ 




— 










"s 


z 










v 


z 










*^ 




600 








--* 


-= 










^"-^_ A 


■■■-.... z 


400 








~~~ ~~ — — -_ __ ___ 


~~ — — — ,5. 


200 


-I 






i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 


■z: 



300 350 400 450 500 550 600 650 700 



Mn 



Figure 1: The energy (in MeV) of the soliton (E so \), 3 times the free-space quark mass 
(3M g ) and the cut-off A plotted as functions of the parameter M (in MeV). Two pa- 
rameterizations of the Gaussian regulator are compared; (i) solid and dashed lines: A is 
fitted to f w = 93 MeV, (ii) lines containing dots: A is fitted to /^ = 1.25 x 93 MeV. For 
the first parameterization E so \ < 3M q at M = 276 MeV and the soliton becomes stable; 
free quarks do not exist beyond M = 309 MeV. For the second parameterization the 
corresponding values of Mq are 345 MeV and 387 MeV. 



mass, M q , is not simply the vacuum value of the chiral field, Mo, but is determined by 
the position of the pole of the quark propagator in the vacuum |J, i.e. it corresponds 
to the solution of k 2 + (r(/c 2 ) 2 Mo + m) 2 \ k2= _ M 2 = 0. The solution for real k 2 exists only 
below a critical value of Mo (see Figure P; above this point no stable free quarks exist. 
However, a stable solution can always be found beyond this point provided the quarks 
dress in a spatially non-uniform background chiral field. 



3 Dependence on the shape of the regulator 

The model (P possesses 4 parameters: the vacuum value of the chiral field M , the 
current quark mass m, the coupling constant G and the cut-off parameter A. The coupling 
constant G is fixed from the stationarity condition in the vacuum, A is fitted to the pion 
decay constant /„., and m to m v . We are left with one free parameter M . 

In this section we analyze how different shapes of the regulator affect the result. We 
compare the Gaussian (g), the monopole ([3]) and a modified version of the instanton ([[]) 
regulator. As we have mentioned it is not possible to continue ([[]) to negative k 2 , which 
is needed in order to obtain the valence contribution. We have instead introduced a form 



which is identical to the instanton regulator for k 2 > while its behavior for k 2 < 
is replaced by a real function which approximately follows the real part of (p]) for small 
negative k 2 ("extended instanton"). 

The Gaussian and the monopole shapes lead to practically identical results. We have 
also tried other shapes and found very similar results. The reason is that several properties 
(including the integral that determines /„.) depend mostly on the behavior of the regulator 
for small value of k 2 . If this behavior is r{k 2 ) «1- ak 2 + . . . then a is almost uniquely 
determined by the value of /„-. Hence, all shapes with this type of behavior lead to very 
similar results. The situation is quite different if dr(k 2 )/dk 2 = 0, e.g. for regulators that 
depend only on k A . We do not find stable solutions for such regulators; the energy of the 
Dirac sea is always higher than the gain due to the lowering of the valence energy. 

The form ([I]) has neither of the above behaviors for k 2 ~ but behaves as 1+^ ^ hi -p- . 
It is therefore interesting to check whether stable solitons can also be obtained for such 
a particular form. We indeed find solitons with very similar properties to those obtained 
using (|) or (|). 

Since f n sets the scale, it is also interesting to study how a higher value for /„. would 
affect the results. 



Regulator 


M 

[MeV] 


A 
[MeV] 


m 

[MeV] 


(ft) 1/8 

[MeV] 


Gval 

[MeV] 


-'-'sea 

[MeV] 


-E'sol 

[MeV] 


/ r 2U/2 

fm 


9a 


Gaussian 
Gaussian 


350 

450 


627 

484 


10.4 
15.9 


-200 

-174 


280 
266 


1715 
1275 


1180 
1261 


1.04 
0.96 


1.16 
1.12 


monopole 
monopole 


350 

450 


834 
639 


5.24 
7.56 


-252 
-223 


275 
260 


2201 
1628 


1176 
1261 


1.05 
0.98 


1.30 

1.28 


extanton 


350 


611 


4.77 


-260 


300 


2374 


1189 


1.05 


1.04 


Gaussian* 


450 


759 


8.75 


-246 


336 


2121 


1458 


0.83 


1.14 



Table 1: Properties of the self-consistent soliton solutions for different shapes of the reg- 
ulator, "extanton" stands for the extended instanton regulator with A = 1/p, Gaussian* 
means a Gaussian regulator with A fitted to /V = 1.25 x 93 MeV, (qq) is the one-flavor 
quark condensate. 



4 Regulators with a 3-momentum cut-off 



The calculation can be made much simpler and the problem of analytic continuation to 
negative k 2 avoided if we take a regulator that does not depend on time. Then the states 
do not depend on v and the integration over v in the above expressions can be carried 
out analytically. All quantities still remain finite and the soliton remains stable against 
collapse. The parameters used are the same as for the 4-momentum regulator. However, 
the explicit breaking of Lorentz invariance leads to serious problems which we discuss in 
this section. 
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Figure 2: a) Different shapes of regulators for M = 350 MeV and A fitted to f n = 93 MeV, 
"extanton" stands for the extended instanton regulator with A = 1/p. b) Comparison of 
the self-consistently determined chiral field for the three regulators displayed in a). 
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The energy of the soliton reduces to 

£soi = N q e val + ]T (e j - ef) + -^ f d 4 x (S 2 + P 2 a 

and the self-consistency equation to 



(13) 




+n c E ?] 







r)< i^rf h 2(( - V2)/A2) ^ (r1 

e,<0 \. l Plh T a,l a ) 



The simplest form of the regulator is a sharp cut-off limiting the 3-momenta to k 2 < A 2 . 
Since one usually uses basis states with good \k\, this is equivalent to restricting the 
Hilbert space. Such a form is very easy to implement since the regulator does not appear 
explicitly in the calculation. Unfortunately, no stable soliton solution is found in this 
case; the contribution of the sea quarks is always larger than the gain due to the lowering 
of the valence energy. 

Stable solutions can be obtained if one takes a smooth form of the regulator, e.g. a 



Gaussian form r(k 2 



-k 2 /2A^ 



. The results are displayed in Figure |3|. The threshold value 



of Mq below which the soliton does not exist is somewhat larger than for the 4-momentum 
Gaussian regulator. One notices a rather striking feature that the soliton does not exist 
beyond a certain value of M . The reason for such behavior is the following: the energy of 



a free quark can be written as e(\k\) = yk 2 + Mlr A (k 2 /A 2 ). For sufficiently small A/M 



the minimum of e(\k\) is not at \k\ = but rather at some \k \ > 0, and the energy of 
a free quark becomes smaller than M Q . Increasing M the cut-off A decreases and at a 
certain value it becomes energetically favorable for the quarks in the soliton to acquire 
sufficiently high momenta and leave the soliton. This happens when 3e(|fco|) < -Ronton- 
Such an unphysical behavior is a clear consequence of breaking the Lorentz invariance in 
the interaction. 

The value of \k \ can be easily determined if we choose a Gaussian regulator and 



m = 0. For the values of A below A c = \/2Mq the minimum of e(| k\) = yk 2 + Moe -2fc2 / A2 
occurs for k 2 = ^A 2 ln(2MQ /A 2 ) and the free quark energy at the minimum is: e( | A; 1 ) = 
A J\(l + ln(2M^/A 2 )). For m/fl this value is slightly modified: 



/c + ln 



2/cM, 2 
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e(N) = A ^ L 2 ^ , /, 
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Figure | shows that the soliton exists for 325 MeV < M < 600 MeV. At M = 
355 MeV the cut-off A reaches the critical value and the lowest free quark state has \k\ > 0. 
Nonetheless, the soliton remains bound since 3e(|/eo|) > Ronton- At M ~ 600 MeV the 
condition is no longer fulfilled and no bound solution exists beyond this point. Between 
these two values the energy of the soliton differs only little from the energy of three free 
quarks. The soliton is weakly bound; the chiral field stays close to its vacuum value while 
the soliton radius is large. 

5 Solitons in the equivalent linear cr-model 

It is well known that the NJL model can be transformed into the form of an equivalent 
linear a-model. This transformation is explained in ||, chapter 5. If one assumes a sharp 
cut-off and sufficiently large A (compared to Mo) the Lagrangian density of the equivalent 
a-model takes the familiar form 

C a = \{d,af + \{d^ a f - ^ (a 2 + n 2 a - flf - \m\ ((a - /.) 2 + ^) . (14) 

The fields a and n a are related to the two components of the chiral field of the bosonized 
NJL model as a(x) = S(x)/g, n a (x) = P a (x)/g, where g = M / f w . The parameter A and 
the mass of the a meson are related to the parameters of the NJL model by0 g = f n /M Q , 
A 2 = 2g 2 and m 2 = 4Mq + m 2 . The question remains whether the above assumptions 
are met in the model described in section 2.0 

We assume that the Lagrangian ( |T4] ) describes the sea quarks while the valence quarks 



are treated separately. This is in the spirit of the approaches used in [1C, 11, ft 
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2 The A 2 here is a factor 2 smaller than the one in Q. 

3 In fact, they are not; the value of Mo /A is close to 1 while the smooth form of the regulator may and 
does generate additional terms in the Lagrangian ( |l4| ) and modifies the values of the parameters. A work 
that will take into account the additional terms is in progress. The conclusions in this section remain 
valid but the results should be considered only as qualitative. 
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Figure 3: Solutions obtained with a 3-momentum cut-off: the energy (in MeV) of the soli- 
ton (bold solid line), 3 times the free quark mass (solid line) and the valence contribution 
to the soliton energy (dashed line) are plotted as functions of Mq (in MeV). The energy 
of the lowest free quark state is equal to Mq + m only below Mq = 355 MeV, above this 
value there exist states with \k\ > which have lower energies. The value \ko\ corresponds 
to the lowest free quark energy. A Gaussian form is used for the regulator, the cut-off A 
(dots) is fitted to /„ = 93 MeV. 



However, in all these approaches the valence orbit is not regularized and - as we shall see 
in this section - this brings a qualitative difference with respect to the situation where 
it is regularized. If we agree that a regulator in the quark-quark or, equivalently, quark- 
chiral field interaction has a well grounded physical origin, the valence orbit should be 
regularized in the same way as the orbits in the Dirac sea. 

The valence orbit is determined as in (|8|). The soliton energy then becomes (a' = 
<r - U) ■ 



2 2 



+ g[af + n a r a ])+E. 



E(a, 7T a ) = 3e val + J d 3 f (± [(Va') 2 + m 2 a a' 2 + (Vvr a ) 2 + 

where the meson self-interaction term (the "Mexican hat" ) is given by 

E s ., = j / d¥ (a' 4 + K 2 ) 2 + AUa* + 2a' 2 vr 2 + 4/^V 2 ) . 
Here we have introduced the source terms which explicitly contain the regulator (see fll2[)): 
f = N q Tes- l qtjf)Pr 2 ((u 2 -V 2 )/A 2 )q Vv0 (f), 
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Figure 4: Solitons in linear cr-models: the solid lines show the energy of the soliton 
(bold), the valence contribution to the energy and of the mass of three free quarks plotted 
as functions of the parameter Mq in the case when the valence state is regularized. A 
4-momentum Gaussian regulator is used; A (see Figure [I]) is fitted to f n = 93 MeV. The 
dashed lines show the respective energies when the valence orbit is not regularized and 
the dotted line represents the corresponding stability line. In this case quasi-stable soliton 
solutions can exist which never happens in the regularized case. 



fa = ^ves' 1 i (*0 ^7 5 r a r 2 ((z, 2 - V 2 )/A 2 )^ (f) 



(15) 



The mean-field solution is obtained by solving the self-consistency equations (the hedgehog 
ansatz a(f) = a(r) and vr a (r) = r a n(r) is assumed): 

<h ' <! - mA a'(r) = f(r) + A 2 \{a\r) + 3/.) a\rf + (a'(r) + f n ) n(rf 




--2~ m * 



7r(r)=#(f)f« + A 2 a'(r) 2 + 2f w a'(r) + n(r) 2 ic(r) . (16) 



In Figures |] and [5] we compare the properties of the soliton when the valence orbit is 
regularized and when it is not. Close to the threshold the two solutions do not differ much 
- as it should be - since here the cut-off A is relatively large. For higher values of Mq (or g) 
the behavior is qualitatively quite different; in particular, the energy of the unregularized 
valence orbit soon becomes negative which never happens for the regularized one. 
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Figure 5: Comparison of the self-consistently determined fields in the NJL model ($ NJL ) 
and in two versions of the linear a-model: with ($ eLSM ) and without ($ LSM - nr ) regulariza- 
tion of the valence state. A Gaussian regulator is used with A fitted to f n = 93 MeV. For 
Mq = 377 MeV, the corresponding A is relatively large (578 MeV) and the three (pairs 
of) curves do not differ considerably; for M Q = 500 MeV, A is much lower (440 MeV) 
and the differences become more important: the fields in the unregularized case acquire 
higher gradients and the soliton shrinks. 



Let us finally mention why the regularization is important also in the linear cr-model 
with only valence quarks. Even if there is no instability as in the regularized linear 
cr-model[], one encounters serious difficulties when going beyond the mean field approxi- 
mation. In the cr-model it is possible to use the Peierls-Yoccoz projection of good linear 
momentum, spin and isospin in order to obtain physical nucleon states. It turns out that 
the energy of the soliton after projection is strongly reduced already when the mean-field 
solution is used for the chiral fields 
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If in addition one allows a variation of the chiral 
field profiles, one obtains a solution with an energy considerably lower than the nucleon 
mass. The energy gain is mostly due to a strongly localized chiral field which lowers the 
valence energy Such a strong localization is not allowed when the regulator is used since 
it puts a physical limit on the gradients of the field (through the source term flT5])). 

In the regularized model both the meson and the quark degrees of freedom are used in describing the 
Dirac sea which leads to appearance of an unphysical pole for the c-propagator p, H . 
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